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Random  pattern  generation  processes  can  provide  a 
rich  class  of  models  for  both  the  synthesis  and  analysis 
of  visual  textures.  This  report  describes  a variety  of 
such  processes  including  point  processes,  random  mosaic 
processes,  and  bombing  processes.  Standard  texture 
analysis  approaches  (e.g.,  cooccurrence  matrices)  are 
related  to  texture  models  based  on  these  random  processes. 
Examples  of  several  textures  synthesized  by  computer  in 
accordance  with  these  models  are  presented. 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESE'*J?CH  fAFSPl 

NOTICE  OF  TKAr,^>aTTAL  TO  DDC  ' 

This  tcohsical  re-iort  L <•- i -n’r ' o-. -.-i  . ^ ^ 

• wua  I]  dis 

rcl9,-.so  lA,/  AFR  (7b). 

13  uaiiiiiited. 

A.  D.  ELOSE 

Technical  Information  Officer  .... 


The  support  of  the  Directorate  of  Mathematical  and  Information 
Sciences,  U.S  . Air  Force  Office  of  Scientific  Research,  under 
Grant  AFOSR-77-3271 , is  gratefully  acknowledged,  as  is  the 
help  of  Mrs.  Shelly  Rowe  in  preparing  this  paper. 

Approved  for  public  release » 
distribution  uniiaited,  * 


1. 


Introduction 


Random  two-dimensional  vector  patterns  are  commonly  en- 
countered in  geology,  metallurgy,  oceanography,  and  in  all 
types  of  photographic  image  processing.  For  convenience,  we 
will  refer  to  all  such  random  vector  fields  as  multicolor 
patterns,  even  when  the  vectors  of  the  process  describe 
features  other  than  color. 

We  are  interesi;ed  in  how  patterns  are  formed  and  what 
par^uneters  are  needed  to  describe  them.  We  will  ask  the 
following  questions:  What  are  the  proportions  of  the  various 
colors?  When  we  move  around  the  pattern,  what  are  the  prob- 
abilities of  the  transitions  from  one  color  to  the  next?  What 
shapes  are  the  pieces  that  comprise  the  pattern,  and  how  are 
these  pieces  combined?  We  would  ''  .to  know  the  first  and 
second  order  statistics  of  our  . < .ns,  e.g.,  histograms, 

variogreuns;  and  the  coocurrences  and  differences  of  colors  at 
the  endpoints  of  a dropped  Buff on  needle. 

One  of  the  principal  applications  of  the  study  of  random 
patterns  is  to  the  analysis  and  synthesis  of  visual  textures. 
First  and  second  order  statistics  are  commonly  used  as 
features  for  texture  classification.  If  textures  can  be 
modeled  as  random  patterns,  it  should  be  possible  to  predict 
the  effectiveness  of  such  features  and  to  design  optimal 
features  for  various  tasks. 

The  purpose  of  this  report  is  to  present  information 
about  a variety  of  spatial  pattern  models,  including  point 
processes,  random  mosaics,  and  bombing  models.  We  describe  a 
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number  of  models  for  spatial  point  patterns  in  Section  2.  The 
"most  random"  of  these,  the  Poisson  point  process,  will  be 
used  in  the  construction  of  most  of  the  succeeding  patterns. 

In  Section  3,  we  will  model  piecewise  contiguous  patterns 
by  random  mosaics.  One  of  these  models,  the  Poisson  Line 
model,  has  interesting  Markovian  properties.  Another,  the 
Occupancy  model,  mimics  natural  cell  growth  processes.  The 
use  of  mosaics  to  depict  patterns  is  of  course  ancient.  The 
Chaldeans  were  skilled  mosaicists  by  2500  BCE.  The  Greeks 
further  developed  the  art  and  were  thought  to  have  used  pattern 
books  for  standard  motifs.  The  pieces  of  a mosaic  were  known 
to  the  Romans  as  tesserae  or  tessellae.  A mosaic  having  a 
simple  geometric  design  was  known  as  opus  tesselatum  — giving 
rise  to  the  current  term  "tessellation". 

We  will  discuss  bombing  models  in  Section  4.  The  bombs 
are  geometric  figures  that  are  dropped  onto  a plane.  These 
figures  in  union  are  the  foreground  of  a pattern,  with  the  un- 
covered portion  of  the  plane  forming  the  background.  Many 
natural  two-color  patterns  are  formed  by  bombing  processes 
or  at  least  can  be  modelled  by  them;  e.g.,  bubble  holes  on  the 
surface  of  cement,  leaves  on  the  ground,  pebbles  on  a beach, 
small  stones  on  the  surface  of  asphalt,  etc. 

Some  statistical  properties  of  our  models  will  be  noted 
in  Section  5.  Finally,  in  Section  6,  we  will  suggest  pro- 
jects for  future  research. 

Realizations  of  all  the  primary  models  covered  in  Sections 


2 through  4 will  be  generated  on  the  computer. 


2. 


Point  Processes 


A two-dimensional  Poisson  process,  with  intensity  A,  can 
be  used  to  describe  a random  distribution  of  points  on  the 
plane.  (The  parauneter  A is  called  the  intensity  or  density 
of  the  process.)  This  process  is  characterized  by  the  property 
that  the  expected  number  of  points  in  a region  of  area  A is  AA, 
irrespective  of  the  orientation  or  shape  of  the  region. 

A two-dimensional  Poisson  point  process  can  be  easily 
simulated  on  the  computer.  A random  number  generator  is  used 
to  generate  A pairs  (xQ,y^,  where  are  the  coordinates 

of  a point  in  the  unit  square  (Fig.  1) . These  points  can  be 
thought  of  as  dropped  onto  the  square,  like  raindrops  falling 
on  a puddle,  or  meteorites  striking  the  moon. 

Not  all  point  distributions  are  random;  a point's  nearest 
neighbors,  e.g.,  may  be  closer  or  farther  than  predicted 
under  a Poisson  model.  Examples  of  such  processes  [1]  are 
briefly  discussed  below,  but  will  not  be  used  to  generate  any 
of  the  textures  discussed  in  this  paper. 

la)  Contagious  model;  Particles  are  attracted  to  their 
neighbors.  An  example  of  such  a process  is  the  spatial  dis- 
tribution of  social  animals. 

lb)  Clustering  model  (center-satellite  process) ; Par- 
ticles are  parents  (or  nuclei)  of  fcunilies  of  children 
(satellites) . The  parents  may  be  dispersed  in  a Poisson 
manner  with  their  children  congregating  about  them.  A good 
ex^u^ple  of  such  a process  is  bushes  in  a field  that  reproduce 
by  sending  off  shoots. 
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is  introduced  by  allowing  the  density  of  points  to  be  a func- 


tion of  location;  i.e.,  X = X(x,y).  For  instance,  A might  be 
dependent  on  local  fertility  or  microclimate. 

2)  Inhibitory  model;  Particles  repel  their  neighbors, 
as  if  they  were  animals  defending  their  territory. 

The  distribution  of  a space-time  point  process  may  be  a 
function  of  a time  varying  intensity  A.  Consider  a forest 
that  starts  with  a random  scattering  of  trees.  As  the  forest 
gets  denser,  the  trees  crowd  each  other  out,  since  they  cannot 
occupy  the  scune  space  at  the  seune  time.  Thus  constraints  are 
imposed  upon  the  location  of  newly  introduced  points.  The  dis- 
tribution of  trees  changes  with  time  from  a Poisson  to  an  in- 
hibitory process. 
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3. 


Cell  Structure  Models 


Many  natural  scenes  are  composed  of  regions  that  appear 
textured;  the  regions  themselves  are  often  composed  of  re- 
latively distinct  patches  of  different  colors.  For  instance, 
when  we  view  a terrestrial  landscape  from  the  air  we  see  a 
mosaic  of  different  plant  species  and  land  uses.  Expensive 
agricultural  machinery  fosters  the  regionalization  of  crops 
into  large  fields.  Urban,  suburban,  desert,  mountainous,  and 
water  regions  all  have  distinct  appearances  from  the  air. 

Such  regions  can  be  modelled  by  random  mosaics  [3].  In  this 
section  we  will  describe  a family  of  pattern  generation  pro- 
cesses that  can  produce  random  mosaics. 

Random  mosaics  are  constructed  in  two  steps: 

1)  Tessellate  a planar  region  into  cells.  We  will  only 
consider  tessellations  composed  of  bounded  convex 
polygons. 

2)  Independently  assign  one  of  m colors  to  each  cell 
according  to  a fixed  set  of  probabilities 

m 

Pi,..., Pm;  J^Pi  = 1.* 


*Other  random  coloring  schemes  are  possible.  Examples: 

(i)  Each  cell  has  a constant  color  throughout,  where  the 
constants  are  Independently  chosen  from  a single  normal  dis- 
tribution. (ii)  All  cell  borders  are  given  a thickness  and 
one  color,  and  all  cell  interiors  are  given  a second  color. 


By  this  process,  we  partition  region  A into  subregions 


m 


iw . . . ,A  ; u A.  = A,  where  A.  is  defined  to  be  the  union 
j=l  3 D 


of  all  cells  of  color  j.  The  partitioning  of  A is  the 
realization  of  a random  process  with  the  following  stationary 
and  transition  probabilities; 

(1)  For  all  s0\.,  Pr(s£A^)  = p^  for  i = l,2,...,m. 

(2)  For  all  (s,s')€A,  with  distance  d = js-s'l  be- 
tween them,  Pr (s ' tA^ 1 s6Aj ) = (d)  = p^(l-W{d;)  + 

6^jW(d)  for  i , j=l , 2 , . . . ,m.  W(d)  is  the  probability 
that  any  two  points  that  are  distance  d apart  are 
both  in  the  same  cell,  and  6^^  is  the  Kronecker  delta. 

Cell  structure  models  form  a family  whose  members  differ 
only  in  the  manner  in  which  the  plane  is  tessellated.  We 
will  describe  some  important  members  of  this  family,  start- 
ing with  the  most  random  members  and  progressing  toward  more 
regular  examples. 

The  first  model  we  will  consider  is  a Gaussian  random 
process . It  is  a degenerate  form  of  a cell  structure  model, 
in  that  each  point  of  an  image  constitutes  an  entire  cell. 

Nondegenerate  cell  structure  models  are  constructed  by 
two  types  of  processes;  (1)  line  processes  and  (2)  cell 
growth  processes.  The  first  type  of  process  uses  straight 
lines  to  partition  a region  into  cells.  The  positions  and 
orientations  of  these  lines  can  range  from  randomly  distri- 
buted (Poisson  line  model)  to  regularly  arranged  (Checker- 


board model) . The  second  type  of  process  uses  a spatial 


3.1  Gaussian  Random  Process 


A real  random  vector  process  {Y^  | s^tA}  is  a Gaussian 

random  process  if  for  every  finite  set  of  points  {s^kA,  the 

corresponding  random  vectors  are  jointly  Gaussian  random 

“i 

vectors.  A (stationary)  Gaussian  random  vector  process  is 

completely  defined  (statistically)  by  its  mean  p = E(Yg  ) and 

“i 

auto-covariance  matrix  Z (b) , where  E (b)  = E[(Y  -p) (Y  -p)^ 

®i+b  ®i 

The  process  reduces  to  one  of  white  noise  if 
E (b)  = E(0)  • A(b) , where  A(b)  = 1 when  b = 0,  and  A(b)  = 0 
otherwise . 
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3 .2  Random  Line  and  Cell  Growth  Processes 

a)  Poisson  Line  Model;  L(  ^ »*•  • • »Pm_i ' • • • ' 

Consider  a system  of  intersecting  lines  in  the  plane 
with  random  positions  and  orientations.  Such  a system  when  de- 
rived by  the  following  Poisson  process  possesses  fundamental 
properties  of  homogeneity  and  isotropy.  A Poisson  process  of 
intensity  t/tt  chooses  points  (9,p)  in  the  infinite  rectangular 
strip  [Ose<iT,  -oo<p<oo].  Each  of  these  points  can  be 
used  to  construct  a line  in  the  plane  of  the  form 
XCOS0  + ysin0  - p = 0,  where  p is  the  distance  to  an  arbitrarily 
chosen  origin.  One  can  use  this  process  to  tessellate  any 
finite  region  A into  cells  (Fig.  2a) . It  can  be  seen  that  there 
are  almost  certainly  four  cells  converging  at  each  vertex. 

These  cells  are  then  colored  in  the  manner  previously  described. 

The  sequence  of  colors  obtained  by  sampling  an  n- 
color  Poisson  line  mosaic  at  equal  intervals  is  an  n-state 
discrete  Markov  chain  with  transition  matrix  P(d)  given 
by 
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Pll(d) 
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P2i(d) 
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...  P„„(d) 
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where  d is  the  sampling  interval.  A matrix  of  the  above  form 
is  a necessary  but  not  sufficient  condition  for  randomness. 

The  additional  requirement  for  the  random  mingling  of  the  color 
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cells  places  a further  restriction  on  P(d);  namely  P.  = 

3a  3 p 

for  all  a,6  ^ j.  A Markov  transition  matrix  having  this 
form  possesses  the  following  properties  [4,  5]: 

(1)  Reversibility  - The  probability  that  color  region 
follows  Aj  in  the  sampling  sequence  is  the  same  as  the 

probability  that  Aj  follows  A^.  This  condition  can  be  ex- 
pressed in  terms  of  the  stationary  and  transition  probabilities 
of  the  process: 

p jP^ j (d)  = p^Pj^(d),  where  (d)  = P^(l-W(d>)  + 6^jW{d) 

and  W(d)  = Q-2xd/i\  ^ 

It  can  be  shown  that  this  condition  reduces  to 

P . . (d)  -P  . . (d)  = a constant  = y for  l^i,  j sn,  i?^j. 

11  3 D 

There  are  many  situations  when  a texture  does, not  possess  this 
property.  For  example,  consider  an  aerial  photograph  of  a 
field  where  the  distribution  of  plant  species  is  controlled 
by  the  prevailing  winds  or  predominant  direction  of  sunlight. 

(2)  Lumpability  - The  chain  resulting  from  any  regroup- 
ing or  renaming  of  the  color  states  is  still  Markov. 

(3)  Specif lability  - A random  ra-color  Poisson  line 
mosaic  is  completely  specified  by  3m  parcimeters.  To  specify 
the  area  covered  by  each  of  the  colors  we  need  to  provide  any 
m-1  independent  terms  of  the  stationary  probability  vector 
(P1P2. . .Pj„) . To  specify  the  transition  matrix  we  need  P^^ 

for  all  i and  j.  However  the  following  restrictions  mean  that 
the  entire  matrix  can  be  specified  by  only  one  parameter, 
namely  the  constant  y: 


(i) 


I Pj^;  1 ^ i,  j i m,  i j 
i:  K 


(ii)  (d) -Pj^^  (d)  = a constant  = y- 


The  constant  y is  proportional  to  the  expected  cell 
2 

width  E{w)  = — . An  additional  2m  parameters  are  needed  for 
the  means  and  covariances  of  the  m colors. 


(4)  The  chain  is  irreducible,  aperiodic,  and  all  its 
states  are  recurrent  states,  since  0 < P^^ (d)  < 1 for  all  i 
and  j.  These  properties  are  defined  as  follows.  The  period  of 
state  a of  a Markov  chain  is  the  greatest  common  divisor  of 
all  integers  n ^ 1 for  which  the  probability  of  returning  to 
state  a,  starting  from  state  a,  in  n steps  is  nonzero.  A 
Markov  chain  is  aperiodic  if  its  period  is  1.  A Markov  chain 
is  irreducible  if  all  its  states  are  accessible  from  each 


other.  A state  of  a Markov  chain  is  recurrent  if  when  we 


start  from  the  state  we  will  eventually  return  to  it. 

Two  or  more  contiguous  cells  of  the  same  color  are 
said  to  form  a patch.  We  would  next  like  to  derive  the  ratio 
of  patch  width  to  cell  width.  Consider  the  distribution  of 
length  of  color  a along  a transect.  It  can  be  shown  [4] 
that  the  sizes  of  neighboring  cells  are  independent  for  a 
Poisson  line  tessellation.  The  points  of  intersection  of  an 
arbitrary  line  with  a Poisson  line  tessellation  constitute  a 
Poisson  line  process  with  intensity  2t/tt.  This  means  that  the 
expected  number  k of  lines  crossing  a transect  of  unit  length 
is  2t/'it.  Thus  is  the  sum  of  j independent  values  of  1/k. 
Let  g()Z^^|j)  denote  the  conditional  pdf  of  given  j.  Then 


A run  of  j cells  of  color  a will  occur  if  the  succeeding  j-1 
cells  are  of  color  a and  a cell  of  another  color  follows  to 
terminate  the  run.  The  probability  that  this  will  happen  is 


p ^”^(l-p  ).  The  conditional  pdf  of  when  j is  allowed  to 

Ot  ot  ^ 

vary  is  therefore 


g{^  ) 


a> 

.1, 

|j)p^~^(i-p  : 

a 

3=1 

-2t  » “> 

e « I 

® j=l 

j-1 


Thus  the  expected  width  of  a patch  of  color  a is  l/d-p^^) 
times  the  expected  cell  width. 

Showing  that  a texture  fits  a Poisson  line  model  re 
quires  the  proof  of  a very  complicated  hypothesis.  We  must 
show  that  i)  the  pattern  is  not  directionally  dependent, 
ii)  cells  of  all  colors  are  randomly  mingled,  and  iii)  the 
pattern  of  each  color  when  taken  against  a background  of  the 
remaining  colors  forms  a 2-color  Poisson  lii  e mosaic. 

b)  Occupancy  Model 

The  occupancy  model  is  defined  as  follows.  A 
Poisson  process  with  intensity  X drops  points  onto  the  plane 
Each  of  these  ooints  spreads  out  to  occupy  a "Dirichlet  cell 


consisting  of  all  points  on  the  plane  that  are  nearer  to  it 
than  to  any  of  the  other  Poisson  points.  These  cells  are  con- 
vex Voronoi  polygons  having,  on  average,  six  sides  (Fig.  2d) . 
(It  can  be  shown  [6]  that  nearest  neighbors  can  be  determined 
in  0(XlogX)  time.)  These  cells  are  then  independently 
assigned  colors  as  usual. 

Two  models  that  are  related  to  the  occupancy  model 
will  be  briefly  noted.  A random  triangular  tessellation  can 
be  constructed  from  Delaunay  triangles  [8]  that  have  as  their 
vertices  the  three  nuclei  that  are  equidistant  from  a vertex 
of  a Voronoi  polygon.  Thus  the  intersections  of  the  borders 
of  Voronoi  polygons  are  the  circumcenters  of  the  Delaunay 
triangles  (see  Fig.  3)  . 

The  Johnson-Mehl  model  [7]  is  often  used  to  describe 
metallurgical  surfaces.  This  model  differs  from  the  occupancy 
model  only  in  that  points  are  dropped  onto  the  plane  as  a 
function  of  time;  i.e.,  X = X(t).  These  points  start  expand- 
ing circularly  as  soon  as  they  hit  the  plane.  A point  on  the 
plane  is  assigned  to  the  cell  whose  expanding  border  first 
reaches  it.  Cells  formed  by  this  process  do  not  have  straight 
line  edges  and  are  not  necessarily  convex  (Fig.  4) . Irregular 
configurations  occur  when  late  arriving  points  fall  near  the 
interface  of  two  large  cells. 

The  occupancy  model  simulates  natural  cell  growth 


processes  on  the  plane.  In  contrast,  true  Poisson  line 
mosaics  are  surely  rare  in  nature.  There  are,  however,  several 
good  reasons  for  choosing  the  Poisson  line  model  as  a standard 


with  which  to  compare  natural  mosaic  patterns: 


1 


I; 


(1)  The  Poisson  line  model  is  mathematically  more 
tractable.  We  do  not  have  an  elementary  foirm  for  W(d)  for  the 
occupancy  model  [9] . 

(2)  The  sizes  of  adjacent  cells  are  independent  for 
the  Poisson  line  model,  but  not  for  the  occupancy  model  [4] . 

(3)  We  would  not  expect  cells  formed  by  a natural 
process  to  have  the  sharp  corners  produced  by  a Poisson  line 
tessellation.  But  a natural  texture  for  which  the  variance 
of  cell  size  closely  resembles  that  of  the  Poisson  line  model 
is  perfectly  plausible  [4]. 
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3.3  Regular  Line  and  Cell  Growth  Processes 


a)  Rotated  Checkerboard  Model;  C (b,Pj^, . . . , M , . . . , 

^1'*  V 

This  is  an  example  of  a cell  structure  model  where 
the  cells  have  a uniform  diameter.  A checkerboard  model  can 
be  formed  by  the  following  procedure.  First  choose  the  origin 
of  an  x-y  coordinate  system  on  the  plane  with  uniform  prob- 
ability density.  Then  tessellate  the  plane  into  square  cells 
of  side  length  b.  Next,  this  "checkerboard"  is  rotated  by  an 
angle  chosen  with  uniform  probability  from  the  interval  (0,2it). 
The  cells  are  now  independently  assigned  one  of  the  m tile 
types  as  before  (Fig.  2h) . The  solution  for  W(d)  for  this 
model  is  discussed  in  the  statistics  literature  as  an  exten- 
sion of  Buffon's  needle  problem; 


W(d)  =1  -4d/iTb  + d 6 b 

1 

= 1-2/n  - (4/iT)cos''^(b/d)  - d^/vh^  + (4/ir)(dVb^-l)^; 

b < d s /2  b 

= 0;  d > /I  b. 


b)  Rotated  Hexagon  Model  [10];  H(X,Pj^, . . . ,P^_j^,yj^, . . . ,y^ 

^l'*-*'^m^ 

This  model  is  analogous  to  the  checkerboard  model, 
except  that  hexagons  are  used  in  place  of  squares.  Another 
way  of  viewing  these  models  is  as  follows;  Consider  a system 
of  particles  on  the  vertices  of  a regular  lattice  (Figs. 5,6a). 
Let  these  particles  be  the  nuclei  of  growing  cells.  Cells  will 


grow  unimpeded  in  a circular  fashion  until  they  reach  the 
tightly  packed  state  shown  in  Figures  5,6b.  At  this  moment 
each  circle  has  four  or  six  points  of  contact  with  its  neigh- 
bors, depending  on  the  nature  of  the  lattice.  As  the  cells  con- 
tinue to  grow  these  points  of  contact  will  be  extended  into 
lines,  and  the  equal  circles  shown  in  Figures  5,6b  will  be 
converted  into  the  equal  hexagons  or  squares  shown  in  Figures 
5,6c.  Notice  the  duality  between  a hexagonal  and  triangular 
tessellation  (Fig.  7). 


f 4 . Bombing  Models 

Random  two-color  patterns  can  be  formed  by  bombing  pro- 
cesses. The  bombs  are  geometric  figures  that  are  dropped  onto 
the  plane.  The  sizes  and  shapes  of  these  figures  are  fixed, 
but  their  positions  and  orientations  are  random.  The  location 
of  a figure  is  determined  by  its  center  point  (i.e., 

center  of  gravity) , and  the  orientation  © of  its  principal 
axis.  By  this  process  a fixed  region  A is  randomly  partitioned 
into  Aj^  and  A^  = A - A^^,  where  Aj^  consists  of  that  part  of  A that 
is  covered  by  the  dropped  figures.  We  shall  refer  to  the 
figures  comprising  A^^  as  isotropically  distributed. 

We  are  assuming  translation  invariance.  Hence  the  number 
of  center  points  falling  on  any  subregion  of  the  plane  depends 
only  upon  the  area  of  the  subregion  — not  on  its  shape  or 
orientation.  The  number  of  center  points  falling  on  any  sub- 
region  A has  a Poisson  distribution  with  mean  XA,  where  X is 
the  expected  number  of  center  points  falling  on  any  unit  area 
of  the  plane. 

To  complete  the  specification  of  our  model,  we  will  color 
regions  Aj^  and  A2  in  a Gaussian  fashion,  with  distributions 
and  N(P2»^2^  respectively.  We  will  now  consider 
three  coverage  theorems  for  bombing  processes. 

Let  K be  a randomly  positioned  convex  figure  with  fixed 
orientation  0.  Let  be  another  convex  region.  Let  a(K^,K,0) 
denote  the  area  in  which  the  center  point  of  K can  be  placed 
so  that  K intersects  K^.  The  area  a(K^,K,0)  has  as  its  bor- 
der the  locus  of  points  where  the  center  point  of  K might  fall 


so  that  the  region  K with  orientation  0 just  touches  the  bor- 
ders of  (Fig.  8). 

The  measure  of  all  center  points  of  the  figure  K where  it 

intersects  K is 
o 

2Tf 

m(KQ,K)  = I a(K^,K,e)  dF© 
o 

where  F(©)  is  the  distribution  function  of  the  orientation 
parameter  0. 

Theorem  I (Dufour)  [11] ; Consider  an  infinite  collection  of 
congruent  figures  K independently,  identically,  and  homogeneously 
distributed  over  the  plane.  The  number  of  figures  K intersect- 
ing another  convex  figure  has  a Poisson  distribution  with 
mean  Xm(KQ,K) . 

Theorem  II  (Dufour)  [11] : Consider  an  infinite  collection  of 

congruent  random  figures  of  area  a amd  perimeter  L distributed 

isotropically,  independently,  and  homogeneously  throughout  the 

plane.  The  number  of  such  figures  intersecting  another  convex 

figure  of  area  a^  and  perimeter  has  a Poisson  distribution 

with  mean  X(aQ  + a + LqL/2it)  . 

Consider  a circular  bombing  process  of  intensity  X.  For 

a circle  of  radius  r,  the  probability  P2  that  a point  chosen 

at  random  on  the  plane  is  isolated  is  equal  to  the  probability 

that  there  is  no  circle  within  a radius  r around  the  point  as  its 

2 

center.  Thus  P2  = exp(-irr  X)  = exp(-Xa),  where  a is  the  area 
of  a circle.  The  proportion  p^  of  the  plane  covered  by 
circles  is  the  probability  that  a point  is  not  isolated: 


* ^“P2  “ 1-exp  (-Xa).  Notice  that  if  we  drop  two  circles  of 
area  l/2a  every  time  we  previously  dropped  a single  circle, 
the  area  covered  will  be  unchanged.  The  fact  that  two  small 
circles  cannot  possibly  combine  to  produce  the  same  shape  as  a 
single  larger  circle  suggests  that  the  area  covered  is  indepen- 
dent of  the  shape  of  the  dropped  figures. 

Theorem  III;  If  an  infinite  collection  of  congruent  convex 
figures,  each  of  area  a,  are  isotropically,  independently,  and 
homogeneously  distributed  throughout  the  plane,  then  the  pro- 
portion Pj^  of  the  plane  covered  by  the  figures  is  1-exp  (-X^). 

Proof;  This  follows  directly  from  Theorem  I.  Let  be  a 

point  on  the  plane.  The  number  N of  randomly  dropped  figures  K 
intersecting  the  point  has  a Poisson  distribution  with  mean 
Xa,  where  a is  the  area  of  K.  Therefore 


f(N)  = (Xa)**  e"^“ 


The  probability  P2  that  a point  is  isolated  is  the  probability 
that  N“0: 

P2  = Pr(N=0)  = f(0)  = e“^“. 

The  probability  pj^  that  a point  is  not  isolated  is  therefore 

, -Xa 
1-e 

Theorem  III  relates  the  number  of  dropped  figures  to  the 
proportion  of  the  plane  covered; 


X = "I  In(l-Pj^)  = - I In 


This  is  a useful  equation  for  estimating  the  number  of  particles 
on  a microscope  slide,  parts  on  a conveyor  belt,  trees  in  a 
field,  etc. 

Examples  of  bombing  processes 

a)  Random  Line  Segment  Process  L{t,£) 

Line  segments  each  of  length  i are  distributed 
isotropically  over  the  plane.  The  orientation  of  a line  seg- 
ment is  specified  by  the  angle  0 between  the  x-axis  and  the 
perpendicular  to  the  line  segment.  The  midpoint 
the  line  segment  specifies  its  position. 


Assume  that  we  are  given  two  line  segments  centered 
at  randomly  selected  points  P and  Q,  with  an  angle  0 between 
them.  They  will  intersect  if  P falls  within  a rhombus  of 
side  length  il  centered  at  Q. 


0 


where  a is  between  0 and  2 it . Thus 


(x-Xq)  COS0  + (y-y^)  sin©  = cosa 


-(x-Xq)  sin©  + (y-y^)  cos©  = sina 


The  equation  for  the  interior  of  the  ellipse  is 


r (x-Xq)  cos©  + (y-y 

L "3^ 


-y^)  sin© 


+ 


-(x-x  ) sin©  + (y-y_)  cos©  < 1 
o o 


The  area  of  the  ellipse  is  a = and  the  perimeter  is 

ff4 

L = 4Cj^E  ^ 2Tt,  using  elliptical  integral  tables  for  E. 

From  Theorem  II,  the  number  N of  randomly  distributed 
ellipses  intersecting  a fixed  ellipse  has  a Poisson  distribu- 
tion with  mean 

4CiE  2 

2XCj^tTrc2+  — ] = TrX[Cj^+C2]  ]. 


The  probability  that  an  ellipse  is  isolated  is 


p^  = Pr(N=0)  = f(0)  = expI-TtX  [Cj^+C2l 

A realization  of  an  ellipsoidal  bombing  process  is  sho%m  in 
Figure  9a. 


1 


c)  Circular  Bombing  Process  [12]  C [ 

A circular  bombing  process  is  a special  case  of  the 
elliptical  bombing  process,  where  Cj^  = C2  = c (Fig.  9b)  . 

Any  two  circles  in  the  plane  will  overlap  if  the  dis- 
tance between  their  centers  is  less  than  their  diameter.  The 
probability  p^  that  a circle  is  isolated  is  equal  to  the  prob- 
ability that  there  is  no  circle  center  within  radius  2r  of 
a center  point  placed  randomly  on  the  plane.  That  is. 


Po  = e 


-4iTr^X 


The  transition  probabilities  for  this  process  are 


given  by  Switzer  [13] : 

P2Pii(d)  = (2P2-1)  + 

where  H(d/r)  = \ - (J^2)  + d s.  2r 

•< 

2 for  d > 2r 


Pl2(d)  = (P2/P2)  ^21^'^^ 

P22(d)  = 

P2i(d)  = 1-P2i(d) 


A special  case  of  the  circular  bombing  process  occurs 
when  each  of  the  circles  covers  only  one  pixel . The  process 
then  becomes  one  of  randomly  thrown  Gaussian  noise  (Fig.  9c) . 

d)  Rectangular  Bombing  Process  R (X ,C2,C2 » y2 ' ^1 ' ^2 ' ^2^ 

Rectangles  of  major  axis  2c2  and  minor  axis  2c2  are 
randomly  dropped  onto  the  plane.  The  orientation  of  a rect- 
angle is  specified  by  the  angle  0 between  the  x axis  and  the 


perpendicular  to  the  major  axis  of  the  rectangle.  The  mid- 
point of  a rectangle  specifies  its  position. 


A realization  of  a rectangular  bombing  process  is  shown  in 
Fig.  9d. 

The  equations for  a rectangle  are 

|(x-XQ)cose  + (y-y^)sinel  = c^^ 

|-(x-XQ)sin©  + (y~yQ)  cos0|  = 

The  interior  of  the  rectangle  is  specified  by 

{l(x-x^)  cose  + (y-yo)  sine|  < i;  [ -(x-x^)  sine  + (y-y^)  cosej 

Cl  Cl  C2  C2 

From  Theorem  II,  the  number  N of  randomly  distributed 

rectangles  intersecting  a fixed  rectangle  has  a Poisson  distri- 

1 2 

bution  with  mean  8X(c,c_  + -(c,+c_)  ).  The  probability  that  a 

12  TT  1 2 

rectangle  is  isolated  is: 

p^  = Pr(:4=0)  = f{0)  = exp(-8X(CiC2  + ^(^1+02)^)). 
Rectangles  of  one  pixel  thickness  can  be  used  to  approximate 
Poisson  line  segments  (Fig.  9f) . 


5. 


Statistics 


The  statistical  properties  of  our  models  can  be  obtained 
directly  from  the  stationary  and  transition  probabilities. 

Note  that  we  have  given  transition  probabilities  only  for  three 
models  — Poisson  Line,  Checkerboard,  and  Circular  Bombing. 

We  will  ass\ime  that  each  region  is  colored  by  a Gaussian 
random  vector  process  having  parameters  and  Let 

f(x,y.,Z..)  = gryi pexp  [ - ^ (x-y ) (x-p ) ] 

a)  Histogram 

Let  H(c)  denote  the  probability  that  a point  on  an 
image  has  color  c,  where  c is  a q-element  vector.  For  example, 
for  visual  color 


H(c) 


Then 


Pr 


“red  content  of  point  s = c^ 
green  content  of  point  s = C2 
_ blue  content  of  point  s = c^ 


H(c) 


m 

I p.f(c,y 
i=l  ^ ^ 


Z.  .) 
11 


b)  Cooccurrence 

Drop  a Buff  on  needle  of  length  d onto  a picture.  The 
endpoints  of  the  needle  are  denoted  by  s and  s’.  Pr(  ) de- 

notes  the  probability  that  s lands  on  color  Cj^  and  s'  lands  on 
V 


color  c 


Then 


1 


+ 


9 


c c 


C, 


I pf  (l-W(d))  f [(/)  , , 1(d)] 


+ I P.W(d)  f[(^^),  1(d)] 


(3)  Difference 

Again  drop  a needle  of  lexigth  d onto  a picture.  Let 
P(A)  denote  the  sum  over  all  colors  c of  the  probability  that 
s'  lands  on  color  c+A  and  s lands  on  color  c.  Then 


1?(A)  = I I 

i 3 


P^(seAj^,s' eAj)P^  [color  (s)  = c,  color  (s')  = 


c+A ( seA^, s ’ eA) ] 


The  equations  under  various  assianptions  are  analogous  to  those 
for  cooccurrence. 


(6)  Variogreun 

Again  drop  a needle  of  length  d onto  a picture. 

Let  V(d)  denote  the  mean  squared  color  difference  at  the  end- 
points of  the  needle.  Then 

V(d)  = E[  (color  (s) -color  (s' ) )^  (color  (s) -color  (s '))  ] 

If  there  is  no  correlation  across  region  boundaries  we  get 


For  the  cell  structure  models,  there  is  no  correlation 
across  cell  boundaries,  and  we  have 


V(d)  = + 2..(0)  + (Ui-yj)]PjPij(d) 


+ I l2E^^(d)]PiW(d) 
i 

+ I t2Z^^(0) ]P^[1-W(d) 1 . 
i 


6. 


Discussion 


We  have  presented  a number  of  models  for  random  spatial 
pattern.  They  not  only  provide  interesting  descriptions  of 
pattern  formation  processes,  but  also  suggest  the  following 
projects  for  future  research. 

Given  a pattern,  we  might  attempt  to  determine  the  best 
fitting  model.  Hopefully,  this  will  offer  some  insight  into 
the  pattern's  structure  and  even  into  the  type  of  natural  pro- 
cess that  created  it.  Preliminary  work  in  this  direction  was 
done  in  an  earlier  report  [19]  . If  we  can  derive  the  para- 
meters of  a model  from  a pattern,  then  we  can  neatly  charac- 
terize the  pattern  as  a realization  of  a particular  random 
process. 

We  have  given  the  sets  of  parameters  required  by  each  of 
our  models.  These  pareuneters  could  be  used  as  features  to 
discriminate  between  different  models  or  different  realizations 
of  the  same  model.  Texture  discrimination  problems  are  funda- 
mental to  the  automated  analysis  of  remote  sensing  data  and 
medical  imagery. 

Our  models  could  also  be  used  to  generate  textures. 

Current  computer  generated  imagery  represents  textured  surfaces 
by  white  noise  or  repeated  patterns.  Random  textures  will  add 
realism  to  computer  animated  films.  Realistic  computer  gen- 
erated images  can  also  be  used  for  visual  perception  experi- 
ments. These  will  offer  other  insights  into  texture  percep- 
tion than  those  obtained  through  the  use  of  random  dot  patterns. 


i 

I 

i 

s 

Much  of  the  research  on  these  models  was  done  by  statis- 
tical ecologists  — Matern,  Pielou,  and  Switzer.  They  were 
interested  in  reconstructing  planar  maps  from  coarsely  sampled 
data.  They  compared  estimated  maps  constructed  by  different 
scunpling  schemes,  under  different  source  errors.  These  same 
techniques  could  be  used  to  reconstruct  a coarsely  sampled 
, image . 

In  an  attempt  to  better  understand  these  models,  para- 
meters that  cannot  be  determined  analytically  might  be 
obtained  by  a Monte  Carlo  simulation.  For  exeunple,  Crain  and 
Miles  [21]  generated  200,000  Poisson  line  polygons  to  estimate 
some  important  unknown  distribution  properties  of  cell  ver- 
tices, perimeter,  and  area.  A similar  study  could  be  under- 
taken to  obtain  the  transition  probabilities  and  cell  statis- 
tics of  the  occupancy  and  Johnson-Mehl  models. 

A number  of  other  models  for  random  spatial  pattern  are 
available.  A broader  study  is  needed  to  compare  and  catalog 
them.  Mandlebrot  [14]  models  irregular  and  fragmentary 
natural  patterns  by  Brownian  surfaces  (Fig.  10).  Wong  [16], 
Hassner  [17] , and  others  model  images  by  two-dimensional  ran- 
dom fields.  Matheron  [18]  and  a number  of  other  French  mathe- 
maticians describe  random  spatial  patterns  in  terms  of  re- 
gionalized random  variables.  Further  study  of  such  models  in 
connection  with  texture  analysis  and  synthesis  would  be  de- 
sirable. 

!! 

S _ _ 


Figure  1.  A realization  of  a two-dimensional 
Poisson  point  process  with  A = 50. 


Figure  2.  Cell  structure  models. 

a)  Poisson  line  tessellation 

b)  L (16,  (.  33, . 33, .34) , (15,30,45)  , 

(25,25,25) ) 

c)  2b  mean  filtered  over  a 3-by-3 
neighborhood 

d)  Occupancy  model  tessellation 

e)  0(40, (.33, .33, . 34), (15,  30, 45), 

(25,25,25) ) 

f)  2e  mean  filtered  over  a 3-by-3 
neighborhood 

g)  Rotated  chec)<:erboard  tessellation 

h)  C(^7,  (.  33  , . 33,  . 34)  , (15,30,50)  , 

(25,25,25) ) 

i)  2h  mean  filtered  over  a 3-by-3 
neighborhood 

j)  Rotated  hexagonal  tessellation 

k)  H (40, (.  33,.  33,.  34) , (15,30,4  5)  , 

(25,25,25))  , o w -j 

l)  2.k  mean  faltered  over  a 3-by-3 

neighborhood 


Figure  3.  Cells  generated  by  occupancy  model. 

Delaunay  triangles  shovm  as  dashed 
lines . 


Realizations  of  bombing  processes. 

a)  E(25,^  , 5 , 20,  10,  50,  10) 

127  127 

b)  C(25,  7,  20,  10,  50,  10) 

127 

c)  C(40,  1,  13,  10,  53,  10) 

127 

d)  R(30,  15,  7,  20,  25,  40,  25) 

127  127 

e)  S(30,  7,  20,  10,  50,  10) 

127 

f)  R(30,  10,  1,  20,  10,  50,  10) 

127  127 
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processes  can  provide  a rich  class  of  models  for  both  the  synthesis 
and  analysis  of  visual  textures.  This  report  describes  a variety  of 
such  processes  including  point  processes,  random  mosaic  processes, 
and  bombing  processes.  Standard  texture  analysis  approaches  (e.g., 
cooccurrence  matrices)  are  related  to  texture  models  based  on  these 
random  processes.  Examples  of  several  textures  synthesized  by  com- 
puter in  accordance  with  these  models  are  presented. v 
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